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ABSTRACT 



This thesis describes the development of a Fortran IV 
computer program for finding shock-induced stresses in a 
three-dimensional piping system. Discretized equations of 
motion are formed by the finite element method. Input 
motions are support translations specified by response 
shock spectra. Maximum responses in individual modes and 
resulting octahedral shearing stresses at selected points 
are found for shock motions in three orthogonal directions 
Extreme stresses, based on the square root of the sum of 
the squares of the modal stresses, are estimated for each 
input direction. An example problem is analyzed to demon- 
strate the use of the program. 
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LIST OP SYMBOLS 



A single underline on a capital letter denotes a 
rectangular matrix, and a single underline on a lower case 
letter denotes a column vector. Superior dots denote time 
derivatives. The symbols used in the computer program are 
described in Appendix E. 

A Element cross-sectional area 

b^ Mode participation factor, mode r 

C A constant 

D Pipe outside diameter 

E Young* s modulus 

Element generalized force vector 
Components of the generalized force vector 
G Shear modulus 

I Second moment of the pipe cross-sectional area 

about a diameter 

Identity matrix 

J Mass moment of inertia per unit length about pipe axis 

K Assembled stiffness matrix 

K Element stiffness matrix 

Z Length of an element 

M Assembled inertia matrix 

M^ Element inertia matrix 

— B12 Element inertia submatrix, bending 1-2 plane 

0 

— B13 Element inertia submatrix, bending 1-3 plane 

0 

Mj^ Element inertia submatrix, longitudinal motion 



7 



”r 

™L 

nip 

N 

E 

Pr 

Q 

E 

r 

s 

ip6 

t 

U 



V 



V 



r 






W 

W 



L 

P 



w 

z 

Y] 



p 

o 

T 



Element inertia submatrix, torsional motion 
Modal mass, mode r 

Mass of the lagging per unit length 
Mass of the pipe per unit length 
Row vector of shape functions 
Vector of principal coordinates 
Principal coordinate, mode r 

Moment of pipe cross-sectional area lying on one side 
of neutral axis 

Displacement vector 

Radius of gyration 

Base displacement 

Kinetic energy of an element 

Pipe wall thickness 

Vector of absolute displacements corresponding to 
unit base displacement 

Components of displacement vector 
Modal matrix 

Spectrum velocity, mode r 

Eigenvector, mode r 

Weight of lagging per unit length 

Weight of pipe- per unit length 

Absolute displacement vector 

One degree-of-freedom system coordinate 

Modified specific weight 

Dimensionless length coordinate 

Density 

Modal stress 

Shearing stress 
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Toct Octahedral shearing stress 

2 

Spectral (diagonal) matrix 

w Natural circular frequency, one degree-of~freedom 

system 

Circular frequency of mode r 
Superscripts 

T Transpose of matrix 

(r) Mode designator for eigenvectors (r=l , 2 , . . . ,n) 
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I. INTRODUCTION 



With the greater demand for electrical power, and the 
diminishing fossil fuel supply, there has been an Increase 
in the construction and use of nuclear power generating 
plants. Correspondingly, there has been greater concern for 
the safety of the power plant during an earthquake. 

Also, the reduction in the size of the Navy has led to 
an increased emphasis on the Integrity of the internal 
systems of a naval vessel under shock loads. Thus, there 
has been a growing engineering interest in finding a rapid 
and accurate method of analysis to determine the shock- 
induced stresses in a complex continuous piping system. 

This thesis describes and presents a computer program, 
written in Fortran IV computer language, to find the 
stresses in a piping system responding to shock. 

Complex piping systems can be analyzed by using dis- 
cretized models. In using the discretization technique, the 
continuous piping system is modeled by finite size elements 
connected at nodes. After the system has been subdivided, 
the elastic and inertial matrices of each element can be 
found and assembled to form the elastic and inertial 
matrices of the system [1,2]. 

A considerable volume of material is available to assist 
in analyzing structures and systems responding to 
earthquakes [3j^,5]. Two of the techniques used in earth- 
quake response analysis are modal analysis with a 
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discretized model and transfer functions [4,5]. The shock 
input for earthquakes can be specified by using time 
history [Sljor shock spectrum [^,5]. 

The Navy uses a technique called the Dynamic Design- 
Analysis Method (DDAM) to determine the shock response of 
shipboard equipment and systems [6]. The DDAM is a modal 
analysis technique using shock spectra to specify the shock 
inputs . 

The shock spectrum presents, as a function of the system 
natural frequency, the maximum response of a one degree-of- 
freedom system responding to a shock motion. The shock 
spectrum is generally presented in terms of response velocity 
or acceleration. 

Several graduates of the Naval Postgraduate School have 
written theses on determining natural frequencies and mode 
shapes for piping systems. Fink [7] developed a program to 
analyze planar piping systems. Including out-of-plane 
bending, using transfer matrices. Baird [8] presented the 
necessary theory to expand Pink*s work to a three-dimensional 
piping system. Rudolf [9] developed a program that deter- 
mines the frequencies of a general three-dimensional piping 
system. Because the transfer matrix does not yield a dis- 
cretized model of the structure, there is no ready means 
for using the frequency and mode shape data to find shock- 
induced modal responses. On the other hand, the finite 
element technique provides a consistent discretized model 
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which may be used for finding frequencies, mode shapes, modal 
responses to shock inputs, and resulting stresses. For this 
reason the finite element method was chosen for this thesis. 

Shock spectra and mode participation factors are used 
to determine the modal responses. At selected points the 
modal octahedral shearing stresses are combined by scalar 
addition of the distortion energy to estimate extreme 
stresses . 
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II. THEORY 



Material and structural linearity and material local 
homogeneity and isotropy are assumed In developing the 
mathematical model. 

A. MODEL LIMITATIONS 

There are few inherent restrictions on the capabilities 
of the finite element method to model complex details of 
piping configurations. Considerations of the quantity of 
input and output data, program length, core storage capac- 
ity, computing time, and cost do impose practical limita- 
tions. The limited time available in developing this 
thesis has necessitated the following restrictions on the 
model used here. 

The system consists of straight lengths of constant- 
section pipe (elements). Each element centerline is parallel 
to one of the three mutually orthogonal global refer- 
ence axes, and 90^ bends are replaced by fictitious exten- 
sions of the tangent sections to the Intersection point of 
the centerlines. Effects of added mass contributed by 
external lagging are represented, but no provision is made 
for added mass due to valves or fittings. Pipe hangers 
furnishing uniaxial restraint in a global direction may 
also be represented. Bending action is described by the 
Euler-Bernoulli beam theory, neglecting shear deformation 






and rotatory Inertia. All ends of the configuration are 
treated as fixed. The shock Input motion Is a translation 
of the base^ along a global axis. 

B. FINITE ELEMENT METHOD 

Once the system has been subdivided Into elements, 
the elastic and inertial characteristics are determined. 

The element stiffness matrix can be determined by many 
different techniques, whereas the inertial matrix is gener- 
ally determined by using shape functions and integrating 
along the length, [1,2]. 

1. The Element Displacement Vector 

Consider an element with displacement components at 
the nodes (ends) as shown in Fig. 1. 




FIG. 1 

ELEMENT GENERALIZED DISPLACEMENTS 

^The clamped ends of the configuration are treated as 
joined to a rigid base. 
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Here, the double arrowhead vectors are positive rotations 
in accordance with the right-hand rule. The displacement 
vector is 



T 

q = [u^ U 2 . . . u^ 2 ^ 

From Fig. 1, components u^ and u.^ are longitudinal dis- 
placements, components u^, u^ , u^, u^, Uq, u^, u^^ and u^^ 
are bending displacements and rotations in the 1-2 and 2-3 
planes, and components u^^ and u^q are twisting rotations. 

2. Element Stiffness and Inertial Matrices 

Przemieniecki [2] forms the element stiffness ma- 
trix by solving the differential equations of the element 
displacements. Appendix A reproduces the element stiffness 
matrix with Przemieniecki ^ s transverse shear deflection 
factor ( 4 >) set equal to zero. 

To determine the element inertia matrix, an expres- 
sion for the element kinetic energy T^ is formed, rotatory 
inertia due to bending being neglected. 

T® = I Cj /q'^ N q d5 = i 4 '^ M® q 

0 - - i - 

where i is the element length, C is a constant/ is the 
element inertia matrix, N is the shape function row vector, 
and c is a dimensionless length coordinate. 
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Therefore 



M® = C£ / N dc 
0 



( 1 ) 



Since the local element axes coincide with the cross 
section principal axes, the inertia matrix may be deter- 
mined from 2x2 and 4x4 element sub-matrices. Consider the 
following motions of the element: 

a) longitudinal: 



and C = pA where A is the cross-sectional area of the 
element and p is the density. From Eq. 1 



and C = J where J is the mass moment of inertia per unit 
length about the pipe axis. Prom Eq. 1 



^ = [u^ = 1 - 5, = C 




I i 



(la) 



b) twist: 





(lb) 



c) bending in the 1-2 plane: 




Ng = t(e-2c^+C^) , Ng = N^2 “ 
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and C = pA. From Eq. 1 




156 22t 54 -13<l 

22l 13JI -3J.^ 

54 13i 156 -22H 

-131 -3l^ -221 4i^ 



(Ic) 



Bending in plane 1-3 results in the same numerical coeffi- 



rotation angles are reversed, rows 2 and ^ and columns 2 
and H are multiplied by minus one. These submatrices are 
assembled to form the element inertia matrix. This matrix 
is displayed in Appendix A. 

C. ASSEMBLED EQUATIONS OF MOTION 

The governing equations of motion, in matrix notation, 
for undamped free vibration [2] with a fixed base are 



where q is the assembled displacement vector and M and K 
are the assembled inertial and stiffness matrices. M and 



are referred to local coordinates. The assembly process 
consists of realigning the local system to the global 
system and building the M and K matrices. Reference [2] 
gives sample assembly techniques. 



dents as • Because the positive senses for the 



M q + K q = 0 



( 2 ) 



e e 

K are referred to a global coordinate system and M and K 
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D. MODAL ANALYSIS 



Modal analysis Is the process of finding the response 
motion of a many degree-of-freedom system by superposition 
of the response motions In the principal (or natural) modes 
This technique is advantageous because the motions In the 
principal modes are Independent of one another (uncoupled). 
To use the method. It Is first necessary to find the nat- 
ural frequencies and mode shapes for the free motions 
governed by Eq. 2, i.e., to solve the eigenvalue problem 
K V = M V. 

2 

Finding eigenvalues (w ) and eigenvectors (v) Is a 
standard problem of numerical analysis. The first step is 
a transformation to a single-matrix problem. Details of 
the transformation used here, which requires a Cholesky 
decomposition of the inertia matrix, are given on p. 3^9 
of reference [1]. Following this, Jacobi rotations [15] 
are used to find the spectral and modal matrices (fi and V) 

The modal matrix is normalized with respect to the inertia 
T 

matrix, i.e., V M V = ^, where is the nth order identity 
matrix. 

The shock analysis is based on using velocity shock 
spectra to characterize the input motion and mode parti- 
cipation factors to determine the responses. Although both 
shock spectra and mode participation factors have been 
extensively used in earlier studies [6], the finite element 
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discretization necessitates a new development of working 
formulas. Details of this development are given in 
Appendix B. 



Each shock input motion consists of base translation 
parallel to a global axis. It is shown in Appendix B that 
the mode participation factor b^ for mode r is 



absolute displacements corresponding to unit base displace- 
ment. The maximum relative displacements due to mode r 
response are 



where is the modal circular frequency and is the 
spectrum velocity at that frequency. 

Since responses to shock inputs in each of the three 
global directions are separately determined, the vector u 
and the mode participation factors b^ must be evaluated 
separately for each direction. The modal masses m^ (see 
Appendix B) are likewise different for each input direction. 

E. GENERALIZED FORCES AND STRESSES 

The subvector q^ of element peak modal displacements 
may be found from the system vector, taking into account the 




(B.5) 



(r ) 

where v is the modal eigenvector and u is the vector of 



q 

—max 




(B.6a) 
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orientation of the local reference axes. From this the 

0 

element generalized force vector f is determined. This 
force vector consists of the elastic and Inertial 
contributions for each element. 

f® = (K® - M®) q® (3) 

Fig. 2 shows the positive directions for the components 
of the generalized force vector. 

f^ = [f f . . . f 

- Li^ I 2 • • • ^12-^ 




FIG. 2 

ELEMENT GENERALIZED FORCES 
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The stresses at the nodal points of the element can be 
determined from the generalized forces by standard methods 
of solid mechanics [10]. 
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III. PROGRAM DEVELOPMENT 



The computer program was written in Fortran IV com- 
puter language using double-precision arithmetic on an 
IBM 360/67 digital computer. It consists of a main calling 
program and seven subroutines. The main program calls the 
subroutines, where the actual calculations occur, in the 
necessary order. Information is passed between the main 
program and the subroutines by the use of a common storage 
core. Appendix E contains the program listing. 

A. SYSTEM DESCRIPTION 

After the piping system has been modeled and subdi- 
vided, subroutine INPUT is used to read the data input of 
the geometry and material properties of the system. The 
geometry of the piping system model is determined by the 
global coordinates of the nodes of the model. The global 
system is a right-handed orthogonal coordinate system. 

Each pipe segment must be parallel to a global axis, but 
there are no restrictions on the placement of the origin of 
the global system. The parameters from the subdivision of 
the system are read. For each node the global coordinates 
and boundary conditions are read. If the node is fixed, 
the stiffness and inertial matrix components due to the 
node displacements are not assembled into the system 
matrices . 
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The numbering of nodes and elements Is arbitrary. 

There are two nodes per element and six degrees of freedom 
per node. For each element, the node numbers and pipe 
group are read. A pipe group includes all elements that 
have the same Young* s modulus, Poisson*s ratio, specific 
weight, radius of gyration, and cross-sectional dimensions. 
Appendix D shows the method used to calculate the effect 
of lagging on the specific weight and the radius of 
gyration. 

B. STIFFNESS AND INERTIAL MATRIX FORMULATION 

Subroutine FGfRM generates the element stiffness and 
inertia matrices as shown in Appendix A. Rotatory inertia 
and shear deformation due to bending are neglected. The 
orientation of the element is determined with respect to 
the global coordinate system and the correspondence between 
local and global degrees of freedom is established. After 
determining the correspondence between the local and global 
degrees of freedom, the element stiffness and inertia 
matrices are assembled to form the system stiffness and 
Inertia matrices. If pipe hangers are present, the corres- 
ponding (uniaxial) hanger stiffnesses are added to the 
appropriate diagonal elements of the system stiffness 
matrix. 

C. FREE VIBRATION MODE SHAPES AND FREQUENCIES 

Subroutine CH£5 mQD uses Cholesky decomposition of the 
system inertia matrix and coordinate transformation to form 



a single symmetric matrix for which eigenvalues and eigen- 
vectors are found. Subroutine JACKET is a modification of 
a program originally coded by Professor G, Cantln. It 
uses the Jacobi variable threshold method to find eigen- 
values and eigenvectors of a real symmetric matrix. The 
modal matrix, when transformed back to system coordinates, 
is normalized with respect to the system inertia matrix. 

Since the system matrices Increase in size with finer 
subdivisions or more complex systems, an economizing tech- 
nique, to reduce the number of degrees of freedom in the 
eigenvalue problem, was investigated during the program 
development. The component mode synthesis method as used 
by Benfleld and Hruda [11] was studied. It was established 
that the constrained component branch technique [11], with 
interface loading, with no node suppression gave results 
identical with those obtained by the direct method described 
above when applied to planar vibration of a clamped-clamped 
beam represented by six elements. Despite the success of 
this trial, it was concluded that the added program com- 
plexity accompanying the use of component mode synthesis 
would offset the potential economies. Another standard 
economizer technique [l4] was considered, but ultimately 
rejected for the same reason. 

D. MODAL RESPONSE 

The number of modes to be used in the stress analysis 
and the specification of the shock spectrum are input data 
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for subroutine SPCTRM. There are three choices available 



to select the number of modes to be used: a designated 
upper frequency limit, an upper limit which is a multiple 
of the fundamental frequency, and an option for any method 
the user desires. Three choices are also available for 
specifying the shock spectrum: a spectrum defined by 
straight-line segments; a constant velocity, then constant 
acceleration shock spectrum; and an option for any method 
the user desires. The shock input is restricted to a 
translation of the base along a global axis and is speci- 
fied by one shock spectrum and three scaling factors for 
the three input directions. Subroutine Mp^DE modifies the 
spectrum velocity for modal mass, if desired. There are 
three choices available: no correction, log of the modal 
weight correction, and an option for any method the user 
desires. This subroutine also calculates the mode 
participation factors. 

E. STRESSES 

Subroutine STRESS calculates stresses at the nodes of 
the elements. The element peak displacements are deter- 
mined and realigned from global degrees of freedom back to 
the local element degrees of freedom. The element gener- 
alized force vector is then calculated. Consider the 
action of the generalized force vector on the element 
shown in Fig. 3» 
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FIG. 3 

STRESS LOCATION POINTS 

At locations ajb,c, and d, the stress vector acting on the 
cross-section may be resolved into normal and shearing 
stresses. Fig. 4 shows the positive sign convention used 
for the normal and shearing stresses. 

The normal stress a and the shearing stress x at each 
point shown in Fig. 3 S-re given by: 

1. at location a) 

o = -f^/A - D/21 

T = D/4l + Q/2It 
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FIG. ^ 

POSITIVE STRESS CONVENTION 

where I is the second moment of cross-sectional area about a 
diameter, D is the outside diameter, Q is the first mo- 
ment, about a diameter, of the portion of the cross-sectional 
area lying on one side of the diameter, and t is the pipe 
wall thickness. 



2. at 


location c) 


0 


= f^/A + D/21 


T 


= D/41 - fg Q/2It 



3 . at location b) 

a = - f^/A - fg D/21 
T = - f|j D/kl + f^ Q/2It 
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4. at location d) 



a = - f^/A - 
T = f^Q D/41 - Q/2It 

These normal and shearing stresses are then used to find 
the octahedral shearing stresses [12], which are given 
by 

1 

^oct " (2a^ + 6 t^)^ (4) 

For each element the above calculations are performed for 
each mode. At each of the four points, the shearing 
stresses for the individual modes are then combined by 
determining the square root of the sum of their squares. 
This process is accomplished for each input direction. 

This is done for each element in turn, starting with the 
first. 

F . PROGRAM OUTPUT 

All input information is echo-checked. The square 
of the mode frequency and the mode shape are printed for 
each mode. The spectrum velocities and mode participation 
factors are also printed. The octahedral shearing stresses 
at points a, b, c, and d (Fig. 3) of each element are 
printed for the three shock input directions. 
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IV. PROGRAM TESTING 



In order to explore the capabilities of the program 
and verify its integrity, a number of plane test configur- 
ations were studied. These are described in the following 
sections. 

A. TEST PROBLEM 1 

The configuration shown in Fig. 5 has three identical 
uniform runs of equal length between the central junction 
and the clamped edges. 




The mode shapes and frequencies for in-plane vibration of 
this system were studied extensively with a developmental 
program. The present program gave identical results for the 
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in-plane modes. The out-of-plane modes exhibited the 
required symmetric or antl-symmetrlc form. 



Additional program tests on this configuration Involved 
different element numbering and node numbering and six 
different orientations of the global axes relative to the 
structure. Mode shapes and frequencies were unaffected by 
these changes. 

B. TEST PROBLEM 2 

This system, shown in Fig. 6, consisted of a straight 
uniform pipe clamped at the ends and represented by four 
elements . 



The frequencies found for the lower modes showed good 
agreement with exact results. Table I shows the comparison 
of the first three bending frequencies of the finite element 
solution with exact results [13]. All the eigenvectors 
showed the required symmetric or anti-symmetric form. 

Bending modes occurred in pairs having equal eigenvalues 
and with the corresponding eigenvectors representing 
deflections in orthogonal planes. 
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FIG. 6 



TEST PIPING SYSTEM 2 
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Stresses found for this system were studied in detail. 
At each node, the stresses are calculated at two points of 
the cross-section. Except at the clamped ends, there are 
two adjacent elements sharing each node, so that stresses 
at these sections are calculated twice. Complete agreement 
was found between these two sets of stress evaluations. 
Also, identical stresses resulted from shocks in the two 
directions perpendicular to the length. 

TABLE I 

COMPARISON OP FIRST THREE BENDING FREQUENCIES OF TEST 



PIPING 


SYSTEM 2 


WITH EXACT 


RESULTS 


Mode 


Exact 

(Hz) 


Finite 

Element 

(Hz) 


Percent 

Difference 


1 


422.32 


422.65 


0.08 


2 


1164.15 


1174.26 


0.9 


3 


2282.20 


2329.64 


2.1 



C. TEST PROBLEM 3 

The equal- legged configuration of Fig. 7 consists of 
uniform pipe throughout. Mode shapes of this system again 
showed the required symmetric or anti-symmetric form. 
Likewise, the shock-induced stresses exhibited the expected 
symmetry properties. 
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FIG. 7 

TEST PIPING SYSTEM 3 
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V. EXAMPLE PROBLEM 



A. PIPING SYSTEM 



The example system analyzed is a model of the main 
steam piping system on a modern naval vessel. The main 
steam piping from the boiler to the rigid cross-connect 
anchor assembly is modeled. The ends of the piping system 
are clamped, and there are two hangers. The system is 
three-dimensional with no branches. The pipe properties 
are : 



outside diameter 
inside diameter 
Poisson’s ratio 
Young’s modulus 
specific weight 
radius of gyration 



7.625 inches 
6.011 Inches 

0.300 

30 x 10 ^ psi 
0.327 Ib/in^ 
^.04 Inches 



The specific weight is a fictitious value which accounts 
for the weight of five Inches of lagging bonded to the pipe. 
Pig. 8 is a schematic of the piping system. The piping 
section PI is 84 inches long, P2 is 264 Inches long, P3 is 
108 inches long, and P4 is 288 inches long. 

B. DATA INPUT AND OUTPUT 

The input data cards are punched in accordance with the 
instructions contained in Appendix E. The echo-check of 
the data is shown in Table II and is arranged as it would 
appear at the end of the program deck. The system 
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FIG. 8 

EXAMPLE PIPING SYSTEM 



eigenvalues and eigenvectors are listed in Appendix C, 

Table VIII. Tables III and IV show the modal velocities 
and the mode participation factors. The stress output is 
given in Table V. 

This problem required a storage capacity of 17^K bytes, 
and 2.67 minutes of computer time. 
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TABLE II DATA INPUT, EXAMPLE PROBLEM 



PROBLEM NUMBER* I 

NUMBER OF PROBLEMS TO BE SOLVED 



TOTAL 

ELEMENTS 



TOTAL 

NODES 



DEG DF FREEDOM 
PER NODE 



DEG OF FREEDOM 
PER ELEMENT 



TOTAL 

GROUPS 



SYSTEM 
DEG OF 



TOTAL 

FREEDQ 



12 


13 




6 


12 1 


NODE 


X COORD 


Y COORD 


Z COORD 


NODAL BOUNDARY 




INCHES 


INCHES 


; INCHES' 


CONDITIONS 


1 


0.0 


0.0 


0.0 


1 


2 


42.00 


0.0 


0.0 


0 


3 


B4.00 


0.3 


0.0 


0 


4 


B4.00 


0.0 


66.00 


0 


5 


B4.00 


0.0 


132.00 


0 


6 


B4.00 


0.0 


19B.03 


0 


7 


B4.00 


0.0 


264.00 


0 


B 


13B.00 


0.0 


264.00 


0 


9 


192.00 


0.0 


264.00 


0 


10 


192.00 


72.00 


264.00 


0 


11 


192.00 


144.00 


264.00 


0 


12 


192.00 


216.00 


264.00 


0 


13 


192.00 


2BB.00 


264.00 


1 


ELEMENT 


LEFT 


RIGHT- 


GROUP 




NUMBER 


NODE 


NODE 


NUMBER 





9 


9 


10 


13 


10 


11 


11 


11 


12 


12 


12 


13 



OUTSIDE 

DIAMETER 

IN. 



INSIDE 

DIAMETER 

IN.' 



POISSONS SPECIFIC 
RATIO WEIGHT 

LB/CU. IN. 



YOUNGS 
MODULUS 
PS I 



SHEAR 
MODULUS 
PS I 



1.154D 07 



PIPE 

GROUP 



RADIUS OF GYRATION 
INCHES 



I 4. 34 

NUMBER OF'HANGERS: 



GLOBAL DOF 



lADOF 



HANGER STIFFNESS 
LB/IN. 



7 

11 



SPECTRUM 

TYPE 



9.B17477D 05 
9.B17477D 05 



FREQUENCY 

TYPE 



BREAK FREQUENCY 
HERTZ 



CONSTANT VELOCITY 
IN 



iNCHES/SEC. 
100.30 



SHOCK INPUT FACTORS 
X Y 



Z 

1 .000 



1.000 1.000 

MODE WEIGHT CORRECTION TYPE 
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TABLE III MODE VELOCITIES EXAMPLE, PROBLEM 

MODE SPECTRUM VELOCITY 

INCHES/SEC. 



1 


100.00 


2 


100.00 


3 


103.30 


4 


100.00 


5 


100.00 


6 


1 33. 30 


7 


100.00 


8 


100.00 


9 


133.00 


10 


88.87 


11 


84.49 


12 


75.10 


13 


71.39 


14 


58.88 


15 


53.46 


16 


47.11 


17 


42.75 


18 


42.19 


19 


32.52 


20 


30.30 


21 


29.92 


22 


29.23 


23 


27.07 


24 


26.38 


25 


25.76 


26 


25.43 


27 


22.07 


28 


21.76 


29 


17.29 


30 


17.04 


31 


16.10 


32 


15.18 


33 


14.34 


34 


14.14 


35 


13.84 


36 


12.24 


37 


12.16 


38 


12.08 


39 


10.67 


40 


10.29 


41 


10.33 


42 


9.65 


43 


9.32 


44 


8.85 


45 


8.35 


46 


8 .31 


47 


8.02 


48 


7.82 


49 


7.53 


53 


7.14 


51 


6.91 


52 


6.61 


53 


6.51 


54 


5.94 


55 


5.80 


56 


5.43 


57 


5.27 


58 


4.74 


59 


4.58 


60 


4.47 


61 


4.30 


62 


3.91 


63 


3.63 


64 


3.43 


65 


2.97 


66 


2.73 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 



IV MODE PARTICIPATION FACTORS, EXAMPLE PROBLEM 



MODE PARTICIPATION FACTORS 
X Y Z 



2.420D 00 


-5.4330-01 


-2.3760-01 


8.881D-01 


1.9360 00 


-8.5790-03 


6.749D-01 


-2.1510-01 


1.3660 00 


1.983D-01 


-8.6390-01 


-2.7390-01 


1.813D-01 


3.6810-02 


4.4490-01 


1.O02D 00 


2.9500-02 


-4.5250-01 


-2.002D-01 


8.8470-01 


-5.9830-02 


1.342D-01 


2 .9090-02 


-4.1320-01 


3.281D-01 


7.3940-02 


-1.4520-01 


-5.372D-01 


2.0870-02 


-2.2340-01 


-1.760D-01 


-2.7650-02 


3.7010-01 


-1.461D-01 


1.7600-02 


-4.3580-01 


-2.436D-02 


-2.0570-02 


-2.5810 00 


2.101D-01 


6.4130-01 


-2.7710-02 


2.120D-01 


3.6170-02 


6.2730-02 


2.414D-01 


-4.2060-01 


-5.3080-02 


-1. 557D-01 


-1.8360-01 


-8.9290-03 


-2.999D-01 


2.7700-01 


-9.8480-02 


-2.523D-03 


-1.5750 00 


2.7830-02 


3.728D-02 


3.9690-01 


4.6230-01 


-7,001D-02 


-7.5210-01 


1.8320-01 


5.680D-02 


-2.8310-03 


3.3900-01 


2.730D-02 


3. 0590-02 


1. 1220-01 


1.372D-02 


-1.4560-01 


-7.4960-02 


2.O96D-01 


3.5930-02 


5.7760-02 


-5,284D-02 


1.7100-01 


-2.9510-03 


3,018D-01 


1 .8490-02 


-2.0460-01 


6, 857D-02 


-4. 1390-01 


3.6120-02 


1.502D-01 


1.9650-02 


2.4800-01 


1.343D-01 


-3.7780-02 


5.7400-02 


1.741D-01 


-5.5900-02 


-1.2860-01 


8.736D-02 


-7.8550-03 


-2.3340-01 


3,2080-01 


-3.3680-02 


-1.1430-01 


2.057D-01 


1.0940-01 


-7.9530-02 


-1.832D-02 


1 .4720-03 


-1.7040-01 


-5.653D-J2 


-2.6610-01 


3.9710-02 


-2.768D-01 


-2. 8490-02 


-1.5850-01 


-1.338D-01 


4.8110-02 


1.0380-01 


5. 282D- 01 


1.8430-02 


-1. 1140-01 


-3.848D-01 


-1.9360-01 


-4.5490-02 


3.217D-01 


-4.5550-01 


1.7640-03 


-9.813D-02 


-2.3620-01 


5.5900-02 


1 .8100-01 


9 .5520-02 


8.5750-02 


-4.9230-01 


-6.4730-03 


-6.5880-02 


-1.5170-01 


4.6710-03 


-9.2750-02 


4.5190-02 


2 .3120-01 


6.8950-03 


-3.7150-02 


8.8880-04 


-4.6700-02 


-7.1100-02 


2.7250-02 


-1.0660-01 


-8.7440-02 


4. 0900-02 


2.8070-02 


-5.4890-02 


-6.5710-02 


-1. 5690-02 


3.3890-02 


4.8980-02 


-3.5650-02 


8.8510-02 


-1.2770-01 


5.8380-02 


-5.3290-02 


-6.4960-02 


-1.2310-01 


5.0540-02 


1 .6870-01 


8.0280-03 


6.4810-04 


-2.6370-02 


-1.6450-02 


-2.0470-02 


2.8140-01 


-7.3650-03 


-2.7480-03 


-5 .7720-02 


-3.5830-03 


1.9600-02 


-2.5910-02 


-1. 7730-02 


-1. 5410-01 


1.0490-02 


8.5290-02 


-1.5480-01 


-1.8900-02 


1.1470-01 


-8. 0590-03 


2.6220-02 


4.9190-02 


1.9480-03 


-1.5460-01 


2.3570-03 


8. 5170-05 


-1.3750-01 


1.0660-04 


-3.4050-02 


-5.9190-03 


-3. 0030-02 


5.4040-03 


-7 .4280-03 


6.0630-03 


2. 5240-01 


4.6170-05 


3.8080-02 



38 



TABLE V OCTAHEDRAL SHEARING STRESSES. EXAMPLE PROBLEM. PSI 



ELEMENT 


1 








9.271550 


03 


6. 508980 


03 


8.408350 


6.674260 


03 


1.687530 


04 


6.482560 


6.609590 


03 


6. 198620 


33 


2.617210 


4.588570 


03 


8.645530 


03 


3.095510 


ELEMENT 


2 








6.609590 


03 


6.198620 


03 


2.617210 


4.588570 


03 


8.645530 


03 


3.095510 


1. 879580 


04 


7.527770 


03 


7.714060 


4.960030 


03 


7.392360 


03 


1.048080 





1 ELEMENT 


7 






03 j 


1.055150 


04 


3.195620 


03 


03 ; 


4.203340 


03 


5.395320 


03 


03 


5.886870 


03 


3.140880 


03 


03 


4.372620 


03 


5.550520 


03 


i 


1 ELEMENT 


8 






03 1 


5.886870 


03 


3. 140880 


03 


03 


4.372620 


03 


5.550520 


03 


03 


2.656300 


03 


3.213450 


03 


03 


7.424390 


03 


1.318900 


04 



ELEMENT 3 
5.764590 03 
I.88813D 04 
4.458580 03 
1.114120 04 

ELEMENT 4 
4.458580 03 
1.114120 04 
4.561880 03 
7.499480 03 

ELEMENT 5 
4.561880 03 
7.499480 03 
4.543860 03 
6.252800 03 

ELEMENT 6 
4.543860 03 
6.252800 03 
3.613900 03 
1.086730 04 



7.613430 03 
■6.471280 03 
5.761720 03 
4.914690 03 



5.761720 03 
4.914680 03 
9.146560 03 
4.30096D 03 



9.146560 03 
.4.300960 03 
7.596550 03 
4.330590 03 



7.596560 03 
4.330590 03 
4.921470 03 
4.381410 03 



1.776200 03 
7.523890 03 
3.033880 03 
2.673220 03 



3.033880 03 
2.673220 03 
3.568610 03 
6.189240 03 



3.568610 03 
6.189240 03 
4.765240 03 : 
7.081600 03 



4.765240 03 
7.081600 03 
2.472250 03 
7.711130 03 



ELEMENT 9 
7.592110 03 
2.469200 03 
5.776440 03 
4.857170 03 

ELEMENT 10 
5.776440 03 
4.857170 03 
7.218960 03 
9.868170 03 

ELEMENT 11 
6.940940 03 
9.868170 03 
1.043000 04 
3.535350 03 

ELEMENT 12 
1. 043000 04 
3.535340 03 
1.769770 04 
6.153500 03 



1.272420 04 
2.811240 03 
1.010420 04 
2.063100 03 



1.010420 04 
2.063100 03 
6.128630 03 
3.488530 03 



6.083360 03 
3.488530 03 
4.379220 03 
2.305370 03 



4.379220 03 
2.305370 03 
1.275160 04 
2.751500 03 



7.845000 03 
2.892360 03 
3.869280 03 
3.820040 03 



3.869280 03 
3.820040 03 
3.567750 03 
4.405200 03 



4.365370 03 
3.508090 03 
3.201110 03 
8.551300 03 



3.201110 03 
8.551300 03 
4.520530 03 
1.570740 04 



4.586790 03 
1.570740 04 
2.335750 03 
9.553940 03 



2.335750 03 
9.553940 03 
6.004770 03 
1.600620 04 



VI. DISCUSSION 



Some of the limitations of the capabilities of the 
program developed above are considered here. 

A. CONFIGURATION LIMITATIONS 

The restrictions that pipe axes must be parallel to a 
global axis and that finite radius bends are not represented 
provide the principal limitations on the configurations that 
can be modeled. There are no inherent limitations on allow- 
able topologipal complexity. 

Certain additional features of real piping systems that 
are not represented in the present modeling Include added 
mass due to fittings and pipe contents, and partial fixity 
or elastic restraint at ends or intermediate points. 

Despite these limitations, it is believed that a signi- 
ficant fraction of current piping systems can be adequately 
modeled using the present program. 

B. SIZE LIMITATIONS 

For the present purpose, the appropriate measure of 
system size is the number of degrees of freedom of the 
system model. This number n, which bears no direct relation 
to physical size, determines the core storage requirements 

and execution time of the program. For approximate estima- 

2 

tion, storage requirements are proportional to n and 
execution time is proportional to n . Using the data from 
the example problem, one can estimate that a 100 



J 



i 




i 



degree-of-freedom system would require about ^OOK bytes of 
core storage and about 9 minutes execution time. 

Increasing the problem size also increases the round-off 
errors. Because double-precision arithmetic (56 bit mantissa) 
is used throughout, observed round-off effects have been 
found negligible (a maximum of 1 unit in the 6th significant 
digit for n = 66) in all applications to date. The very 
small errors detected are attributed to the residual eigen- 
vector errors present when the Jacobi rotations are 
terminated . 

In view of the foregoing considerations, it is believed 
that the practical upper limit on problem size (with the 
IBM 360 / 67 ) is about n = 110 and is determined by core 
storage capacity. 
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VII. CONCLUSIONS 



It is concluded that an effective program has been 
developed for determining shock-induced stresses in piping 
systems. The principal limitations of the present version 
can be removed by adding features that are clearly within 
the current state-of-the-art. Recommended additions are 
listed below; 

1. Remove the pipe axes orientation restriction so that 
a general piping system can be analyzed. 

2. Develop . element stiffness and inertia matrices for 
bends. 

3. Modify the program so that partial fixity and elastic 
restraint at the ends or intermediate points may be 
included in the boundary conditions. 

4. Include the mass effects of valves, fittings, and 
pipe contents. 

5. Replace the Jacobi rotation method by a more efficient 
eigenvalue-eigenvector algorithm. 
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APPENDIX A 



ELEMENT STIFFNESS AND INERTIA MATRICES 
TABLE VI Element Stiffness Matrix 



EA 

Z 

0 

0 

0 

0 

0 

” z 

0 

0 

0 

0 

0 



12EI 



z-^ 

0 

0 

0 

6EI 
£2 

0 

12EI 



12EI 

£3 
0 

-6EI 

£2 

0 



.3 



6EI 

O 



0 

12EI 
■ .3 

0 

6EI 



2GI 

Z 

0 

0 

0 

0 



4ei 

z 

0 

0 



4EI 

z 



0 

2GI 

Z 

0 

0 



0 -> 

6ei 

„2 



2EI 

Z 



Symmetric 



6E1 

2 
X 

0 

0 

0 

2EI 



EA 

Z 

0 

0 

0 

0 

0 



12EI 

Z^ 

0 



0 

6EI 



12EI 

.3 

0 

6EI 

„2 



I 



i 

J 



f 




I* 



TABLE VII Element Inertia Matrix 



t o* 
: c\j 



: O 

; c\j 




< on 



: o 
: c\j 



; o 



on ' 

iH 



< OJ 
cx cr 



: o 
: c\j 




n 



SI 






APPENDIX B: VELOCITY SHOCK SPECTRA AND MODAL RESPONSE ^ 

A shock spectrum exhibits the maximum response dis- 
placement of a single degree-of- freedom system whose base 
is subjected to the shock motion. Consider a single degree- 
of-freedom whose base has a shock displacement s, and whose 
displacement relative to the base is z. The equation of 
motion is 



z 



2 

+ Z 



s 



(B.l) 



where w is the natural circular frequency of free vibration 

with the base fixed (s = 0). If z represents the ex- 
max 

treme value of z in response to the shock motion s^ then the 
displacement shock spectrum for that shock motion is a plot 
of z versus w. Equivalent information can be given in 
the form of a velocity shock spectrum. In this form, the 
spectrum velocity V is related to the maximum displacement 



V = 



(B.2) 



To develop the equations for finding modal response 
from shock spectrum data, let w be the vector of absolute 



^Thls Appendix is based on material presented by Prof. 
R.E. Newton in the course ME^5223 September 1971- 



^5 



displacements in the system degrees-of-freedom. For a 
uniaxial translation s of the base, this may be expressed 
as 



w = £ + us 



(B.3) 



where q is the vector of displacements relative to the 
base, and u is a vector of constants. To illustrate the 
meanings of these vectors, consider the beam element of 
Fig. 9. 




\\\ Ineri\al Re-fc rcnce \ 



FIG. 9 

BEAM ELEMENT COORDINATES 

T T 

For this example, q = [q^ q^ q^ q^] , w = Cq^+s q^ q^+s q^^] , 

and u = [1 0 1 0]. It can be seen that u represents the 
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I 



ft 












absolute displacements resulting from unit base 
translation. 

The equations of motion of an n degree-of-freedom 
system with base motion may be written 





M w + K £ = 0 


Using Eq. 


B.3a this becomes 
M£+K£=-Mus 


Using the 


substitution q = V £ where £ is the vector of 



principal coordinates, this may be rewritten as 
v'^’mV£ + v'^KV2. = -v'’^ M u s 



or 


£ + g_ = - b s (B.it) 


where b = 


T 

V Mu. Eq. B.4 is equivalent to the scalar 



equations 

Pjo + Pp = - s , (r=l,2, . . . ,n) (B.^ta) 

where p^ and b^ are the rth components of £ and b, 
respectively, and is the circular frequency of mode r. 
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The component is called the mode participation factor 
for mode r. It is given by 




where is the eigenvector for mode r (rth column of 

V). 

Comparing Eq. B.4a with Eq. B.l, and using Eq. B.2, 
the maximum response in mode r may be expressed as 



(p ) = b V /o3^ 

^^r^max r r r 



(B.6) 



where is the spectrum velocity at frequency Using 

Eq. B.6, the corresponding extremum for the relative 
displacement vector is 






^max = ^ 



b V /o3 
r r r 



(B.6a) 



Eq. B.6a is used for calculating peak modal responses in 
order to find stresses. 

In shipboard shock studies, it has been found that 
massive equipment may partially suppress the base motion 
[6]. This may be taken into account by appropriate modifi- 
cation of the spectrum velocity V^. Such a correction re- 
quires apportionment of the system mass among the modes. 

The modal mass for the rth mode is taken to be 
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APPENDIX C: EXAMPLE PROBLEM EIGENVALUES AND EIGENVECTORS 



Table VIII shows the first six eigenvalues and the 
corresponding eigenvectors of the example problem. The 
eigenvalue (square of the circular frequency) appears at 
the head of a column followed by the components of the 
eigenvector grouped by nodes. The remaining eigenvalues 
and eigenvectors are omitted. Spacing of the eigenvalues 

remains approximately uniform throughout the remainder of 

0 

the set. The largest eigenvalue is 1.32 x 10°. 



l\9 



TABLE VIII 



EXAMPLE PROBLEM EIGENVALUES AND EIGENVECTORS 
3.20210D 02 1.76822D 03 4.2839AD 03 7.28442D 03 



1.36643D 04 1.93694D 04 



3.734190-05 
2. 419560-03 
4. 498290-03 
2.613170-04 
-9. 722550-05 
1.227360-04 



1.346630-05 

1.090030-02 

1.041630-03 

-1.297640-03 

-2.247990-05 

3.964450-04 



2.084720-04 
-2.656940-02 
1. 106640-02 
4.111110-04 
-2.557670-04 
-1. 134310-03 



-1. 140300^04 
-1.284070-01 
-5.213800-03 
1.334410-03 
1.123570-04 
-5.432470-03 



-4.651040-04 

8.533350-02 

-1.597670-02 

1.670380-04 

3.099150-04 

3.613510-03 



9.202400-04 

5.243660-02 

2.803910-02 

3.876080-04 

-5.483340-04 

2.206790-03 



7.4 683 00-05 
1.094290-02 
-1.656710-03 
5. 226180-04 
5. 072280-04 
2.907080- 04 



2.693140-05 
2.303660-02 
-3. 865930-04 
-2. 594880-03 
1.175200-04 
6.000790-05 



4. 169000-04 
-8.449300-02 
-1.211330-03 
8. 219070-04 
1.109210-03 
-1.500320-03 



-2.280180-04 

-4.007280-01 

1.910950-03 

2.667080-03 

-5.856580-04 

-6.912510-03 



-9.298920-04 
2.678550-01 
1.144910-02 
3.336680-04 
-2. 056210-03 
4.702150-03 



1.839590-03 

1.628400-01 

-1.906140-02 

7.738690-04 

3.551720-03 

2.824790-03 



7. 586450-02 
-3.272250-02 
-1. 502070-03 
7. 817170-04 
1.680480-03 
6.803640-04 



1.665440-02 
2.268140-01 
-3.508370-04 
-3.299150-03 
3.475920-04 
-1 .369440-03 



1.461590-01 
-1.378830-01 
-8.567620-04 
5.612880-04 
2. 712830-03 
-1. 504380-03 



-7.412730-02 

-5.514730-01 

1.733630-03 

9.105430-04 

-1.334310-03 

-6.234990-03 



-2.396660-01 

2.053020-01 

1.086020-02 

1.877200-03 

-3.910450-03 

4.765840-03 



3.935970-01 
8.261260-02 
-1.804080-02 
1.713970-03 
5.854290-03 
2. 827700-03 



2.077980-01 
-8. 986820- 02 
-1.347380-03 
9.335880-04 
2.214390- 03 
1.069920-03 



4.113490-02 
4.263500-01 
-3. 150070-04 
-2.558230-03 
3. 624590-04 
-2.7978 10-03 



2.873100-01 
-1.358610-01 
-5. 017450-04 
-7.254990-04 
1. 113250-03 
-1. 505580-03 



-1.404710-01 

-4.479070-01 

1.554770-03 

-4.148180-03 

-4.794290-04 

-5.537350-03 



-3.879600-01 

2.414480-02 

1.025300-02 

3.312280-03 

-4.396840-05 

4.800720-03 



5.383140-01 
-4.128320-02 
-1.697730-02 
1. 745820-03 
-2.271960-03 
2.806380-03 



3. 554280-01 
-1. 541720-01 
-1. 192630- 03 
1.005280-03 
2. 171520-03 
1.459330-03 



6. 127200-02 
5.499210-01 
-2. 791070-04 
-1.168140-03 
2.325480-04 
-4.223990-03 



2. 505250-01 
-3.754620-02 
-1.464640-04 
-2. 229600-03 
-2.408720-03 
-1. 503930-03 



-1.256530-01 
-4.379060-02 
1.374510-03 
-7. 483740-03 
8.999310-04 
-4.821860-03 



-2.500440-01 

-1.789320-01 

9.628570-03 

2.470720-03 

3.612420-03 

4.806570-03 



1 .421630-01 
-1.135150-01 
-1. 587330-02 
3.151950-04 
-8.321710-03 
2.761020-03 



4. 844370-01 
-2.219640-01 
-1.037840-03 
1.050640-03 
1.676140-03 
1. 848530-03 



7. 136440-02 
5.922980-01 
-2.431470-04 
-3. 187010-04 
8.266200-05 
-5. 646860-03 



-3. 727910-02 
1. 513220-01 
2.088940-04 
-3.438590-03 
-6.263850-03 
-1.499420-03 



-3.999870-02 

4.44113D-01 

1.193020-03 

-6.623490-03 

1.465560-03 

-4.090820-03 



-5.137650-02 
-2. 818790-01 
8.987980-03 
8.151150-04 
1.029980-03 
4.783360-03 



-3.036070-01 

-9.597510-02 

-1.473150-02 

-5.005620-04 

-3.163490-03 

2.692010-03 



4. 544510-01 
-1. 159680-01 
-7.857740-02 
1. 111580-03 
1.227640-03 
2.067390-03 



7.134900-02 

2.764270-01 

-2.704960-03 

-1.873420-04 

1.566670-05 

-5.781560-03 



-3. 730900-02 
7. 101190-02 
3.940520-01 
-4.581830-03 
-7. 87601D-03 
-1.448660-03 



-3.986070-02 
2.253770-01 
-7.170640-02 . 
-3.766930-03 
1.203200-03 
-4. 120900-03 



-5.067870-02 
-5 .032210-02 
7. 375440-02 
3.146250-06 
-2.844290-03 
3.231510-03 



-3. 040550-01 
2.482560-02 
-3.788760-02 
-7.118030-06 
3.160610-03 
1.327770-03 



4.844510-01 
-5.323830-05 
-1.372120-01 
1. 172420-03 
9. 790640-04 
2.211200-03 



7.132310-02 
1.192160-03 
-2.735590-03 
-5.588410-05 
-6. 817080-06 
-4. 079800-03 



-3.732580-02 
2.462120-04 
8. 116120-01 
-5.719250-03 
-7.355270-03 
-1. 107360-03 



-3.969890-02 

4.440830-04 

-1.282270-01 

-9.022450-04 

9.235360-04 

-4.114530-03 



-4.992390-02 
-2.334730-03 
2.579220-01 
-8. 088350-04 
-3.532320-03 
-2.122720-03 



-3.040190-01 

-1.372800-03 

-2.657550-01 

4.863670-04 

4.543460-03 

-2.715950-03 



3.223770-01 
-3.993090-05 
-5. 555850-02 
1.023740-03 
7.344430-04 
2. 254460-03 



2.421760-01 
8.943870-04 
-2.962940-03 
3.440720-05 
-5.118530-06 
-7. 879640-04 



1.778870-02 
1.847940-04 
3.684660-01 
-6. 240510-03 
-5.531030-03 
-4.405480-04 



2.072240-01 
3.334770-04 
-1.034650-0 1 
1.293160-03 
6.957710-04 
-2.320910-03 



3.027290-01 
-1. 755140-03 
1.527010-01 
-2.121270-03 
-2.671720-03 
-5.611510-03 



6.348150-02 
-1.033020-03 
-1. 749090-01 
2.211040-03 
3.448750-03 
-6 .110440-03 



1.674890-01 
-2. 662160-05 
-6. 529060-04 
4. 239510-04 
4.896970-04 
1. 980800-03 



2.105480-01 
5.963870-04 
-8 .495620-05 
2. 892450-05 
-3.415200-06 
1.424270-03 



2.932840-02 

1.232610-04 

5.640670-03 

-3.012020-03 

-3.694300-03 

8.409130-05 



2.586930-01 

2.225160-04 

-2.713260-03 

9.958980-04 

4.653360-04 

8.950440-04 



5.251100-01 
-1. 172050-03 
3.420560-03 
-1 .428070-03 
-1.791900-03 
3.088030-04 



3. 673960-01 
-6.903080-04 
-4. 391680-03 
1.732520-03 
2.318890-03 
-1.287430-03 



4. 818110-02 
-1. 331100-05 
7.315500-03 
-9.921070-05 
2.448690-04 
1.249900-03 
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1.165380-03 

4.226120-03 
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APPENDIX D: CALCULATION OF MODIFIED SPECIFIC WEIGHT 



AND RADIUS OF GYRATION 

There is negligible effect on the stiffness of the 
pipe element from the lagging on the pipe; however, the 
lagging does contribute significantly to the mass of the 
element. By appropriately modifying the pipe element 
specific weight and radius of gyration, the mass effect of 
the lagging can be included. The modified specific weight 
is 



+ Wl)/A (C.l) 

where is the modified specific weight of the pipe ele- 
ment, and the sum (W + W,. ) is the combined weight of the 

p 

lagging and pipe per unit length. The radius of gyration 
is given by the relation 

J + = (m„ + m, ) (C.2) 

p L p L 

Where J and Jr are the mass moment of inertia per unit 
p L 

length of the pipe and lagging, respectively, (m + m,. ) is 

P L 

the combined mass of the pipe and lagging per unit length, 
and r is the radius of gyration. 
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WRITE(6 ,195) 

DO 160 J=1,IGRPT 

READ( 5, 183) DO ( J ) , D I ( J ) , POI ( J ) , SP WT ( J ) , E ( J ) 
READ(5,197) EFRAD(J) 
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G{ J)=E( J)/ ( 2 .D 0 ^( 1 .DO+POI ( J ) ) ) 

160 WRITE( 6 , 196 ) 0 0 ( J ) t DI ( J ) , P 01 ( J ) t S P WT ( J ) t E ( J ) , G ( J ) 
WRITE( 6 , 198 ) 

WRITE ( 6 , 199 ) ( J, EFRAD( J) ,J= 1 , IGRPT) 
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RETURN 

END 
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DO 3160 I=1,ND0FPE 
DUM( I )=SM( 1,1) 
DUMKI ) =SM(I ,2) 
3160 DOME ( I )=SM( 1,3) 
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